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ABSTRACT 


One of the non-linearities which occur in feedback control systems 
is transportation lag. All of the currently popular mathematical tech- 
niques available for the analysis of these systems are severely compli- 
cated by the introduction of transportation lag. The suitability of a 
numerical method employing the Z-transformation to this problem was in- 
vestigated, as well as the particular advantages and disadvantages of 
the other methods. Solutions by the numerical method were compared 
with solutions by other methods and the results obtained presented in 
tabular form. These results demonstrate the superiority of the numeri- 
cal method in determining the transient response of these systems with- 
in reasonable engineering accuracy in a much shorter time. 

The writers wish to express their appreciation to Dr. M. P. Pastel 
for his guidance during the initial stages of the work, and to Dr. G. J. 
Thaler for his advice and suggestions throughout the investigation and 


his aid-in the preparation of this paper. 
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aie Introduction 
1.1 Purpose of the Study 

Feedback control systems are seldom truly linear, but rather 
they are subject to the introduction of a considerable number of non- 
linearities. Some non-linearities may be introduced intentionally by 
the system designer. An example of this might be the selection of a 
relay servomechanism with its attendant non-linearity, when a linear 
system could accomplish the task just as well. On the other hand there 
may be no way of avoiding a non-linearity in the system design. An ex- 
ample of this would be the mandatory inclusion of a gear train with its 
attendant backlash. 

A non-linearity which frequently can neither be avoided nor 
PST ey is transportation lag. 

Because of the variable influence of transportation lag in 
feedback control systems and its avoidance whenever possible, it has 
received relatively little prominence in the literature. Rather, the 
mathematical tools capable of analyzing systems with transportation lag 
have been indicated without further amplification. While it is true chat 
the tools indicated for the analysis of systems with transportation lag 
will yield satisfactory results, the problem does not end there. These 
standard tools, while yielding accurate results with a reasonable amount 
of work for linear systems, require considerably more work and time when 
the problem is complicated by transportation lag. Consequently, a need 
for a better tool than those currently in use for the solution of pro- 
blems involving transportation lag is indicated. 


In 1954, a paper by Ragazzini and Bergen outlined a numerical 
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methods technique for determining the transient response of linear feed- 
back control systems. The nature of the mathematics which the method 
used indicated that it might handle the non-linearity of transportation 
lag more conveniently than the aforementioned standard mathematical 
methods generally known to engineers. The purpose of this study is to 
evaluate the numerical methods technique of Ragazzini and Bergen in 
analyzing feedback control systems with transportation lag and to com- 
pare the accurecy attainable and the work involved in using this method 
with the accuracy attainable and work involved in using more popular 
methods. 

In doing so, it will be necessary to describe briefly the 
common analytical and graphical methods available for solving linear feed- 
back systems. The increased difficulty and complexity introduced by the 
transportation lag into these methods will be demonstrated. Finally, 
the application of the numerical method to the same problem will be 
presented. 

1.2 Summary of the Results and Conclusions 

The numerical methods technique was compared with the root 
locus, Bode diagram, power series expansion, and second order Pade approxi- 
mation approaches to the problem of system analysis. The Bode diagram 
proved to be the speediest method of all for determining system stability, 
but this method yields no other specific information regarding the over- 
all system performance. For the full transient response the numerical 
methods technique proved to be faster and more accurate than the second 
order Pade approximation. The root locus method, while retaining its 
potential for good accuracy, required a prohibitive amount of time and 
tedium to be considered a practical working tool for problems involving 
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transportation lag. The power series expansion approach to the pro- 
blem proved to be totally impractical for other than extremely simple 
systems. 

In short, the numerical methods technique is the fastest and 
most accurate means of determining the transient response of a feedback 
control system with transportation lag. Moreover, the method can be 
understood and applied by engineers and students whose background includes 
only a first course in linear servomechanisms. In the solution of the 
problem by the numerical methods technique, no more complicated computing 
device than the desk calculator is required. 

In all cases, solutions determined by the various methods were 
compared with the transient response obtained from programming the Fourth 
Order Pade Approximation of the Laplace shift operator on a ten channel 
Donner analog computer. This result was the most correct obtainable 


under the conditions of time and equipment available for the study. 
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Ze Description of the System 


The feedback control system to be studied is shown in block 
form in figure 2-1. In order to carry out the goals of this thesis, a 


number of requirements were set for this system. 





A frequency response diagram crossover of one radian per 





second for all variations of the system was desired. For this reason, 


| 


| the value of K in block one will always be equal to one. 

The forward loop transfer function must always be low pass in 
nature. For this reason the following values of G(s) were selected in 
, block two. These values also give us convenient first, second, and 


third order systems to investigate. 


G(s) = 1 Zed 
| 
G(s) = ——-- 222 
o O.1S +/ 
G(s) = ! eG 


(0.18 +1) (0.055 +1) 

Block three contains the transportation lag. The nature of 
this block and its analytic representation will be discussed in detail in 
section 4. 


The system will use unity feedback. 
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SB Linear System Analysis 


3.1 Discussion 

Various methods of analysis and solution of linear feedback 
Systems are available. Two graphical methods currently in use supply 
the necessary information required for the design and analysis of these 
systems. For the purposes of this study, the transient, or time varying 
response characteristics, are more important. 

3.2 Frequency Response 

The Bode network therorems presented a new approach to the 
Standard Nyquist stability criterion. One of the great values of the 
Bode diagram is the simplicity with which approximate results of suf- 
ficient accuracy may be obtained. 

The Bode diagram for the linear second order system described 
in section 2 is shown in figure 3-1. Phase margin and gain margin values 
are indicated, and reasonable predictions of system performance may be 
made from these values. This graphical presentation of the open loop 
frequency response is plotted knowing the location of the open loop poles 
and zeros. 

System stability may be inferred from the values of gain margin 
and phase margin on the open loop frequency response. Information from 
the Bode diagram may be transferred to the Nichols chart, and the closed 
loop frequency response may be obtained. Information from the Nichols 
chart should agree closely with that from a basic Nyquist plot. 

3.3 Root Locus 

A graphical method exists for determining the locus of the 

roots of the closed loop system from the location of the poles and zeros 


of the open loop system. The root locus method essentially demonstrates 
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Frequency Response Diapram of Undelayed Second Order System 


Figure 3-1, 
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the transient response characteristics of the system rather than the 
frequency response. 

The root locus plot for the linear second order system describ- 
ed in section 2 is shown in figure 3-2. Since this method indicates 
directly the effect on the closed loop response of the system gain from 
zero to infinity, it is obvious that this system is stable for all values 
of K. This may also be inferred from inspection of figure 3-1, which 
shows that the phase angle curve always remains above the 180 degree 
phase angle line. 

3.4 Transient Response 

The transient response of a control system contains nearly all 
the parameters required by the engineer for the evaluation of its perfor- 
mance, A typical transient response is shown in figure 3-3 and the vital 
parameters are indicated. From a study of these values, an estimate of 
the worth of the control system may be made, and the need for compensa- 
tion may be determined. 

The parameters which will be of interest are the following: 
Rise Time - Length of time for the system to first reach a point within 

5% of the desired value. 
Peak Overshoot - Value in percent by which the system overshoots the de- 
sired value. 
Time to Peak Overshoot - Length of time to reach peak overshoot point. 
Settling Time - Length of time for the system to settle within 5% of the 


desired value. 
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Figure 3-2. Non-Dimensionalized Root Locus of Linear Second Order System 
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tee Transportation Lag 


4.1 Discussion 

"Transportation lag'' is a term frequently used in the process 
industry to describe a pure time delay. Other terms used to describe 
this phenomenon include ''dead time", "delay time", and "true time lag". 
Figure 4-1 outlines an arrangement which might be used to control the 
thickness of steel plate produced by a rolling mill. A sensor indicates 
the actual plate thickness in voltage form, The difference between its 
output and another voltage representing the desired thickness forms the 
error voltage. The amplified error ultimately positions the rolls, The 
transfer function from plate thickness at the sensor X through the ampli- 
fier, motor, and gear train to roll position, hence thickness at the rolls 
Y, differs in no sense from ordinary transfer functions considered pre- 
viously. From Y to X however a new element is introduced because a finite 
time must elapse before a change in thickness at Y reaches the point of 
measurement at X, d units of length away. If the steel plate proceeds to 
the left with a constant velocity v, any time function of thickness y(t) 
resulting from changes in roll position will be reproduced at the sensor 
without distortion d/v time units later. It is precisely this time quan- 
tum that the term transportation lag describes. 

Transportation lag may be expressed by the following mathe- 
matical relationship: Given f(t) for t =0, the function f(t-T) differs 
from the function f(t) by the dead time T, where f(t-T) is defined to be 
zero for t<T. 

The setup of problems on computers occasionally requires the 
representation of transportation lag or an expression of the form f(t-T). 


The brute-force approach to the solution of this problem, and one that 
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has been used in the past, is to generate f(t) in the computer and plot 
f(t) vs. t on a specially equipped recording device. Schemes involving 
such mechanical devices to accomplish the introduction of transportation 
lag have proved to be workable for some problems, but the limitations 
imposed by the physical design of the devices and the introduction of 
human error in their use have made necessary the development of better 
methods, 


The Laplace transform of f(t-T) is 


LCst-T)] = Lf sees” = Fares’ - 


if f(t-T) = 0 for t<T. In equation 4.1 4 is commonly called the 
Laplace shift operator. The problem of generating a function f(t-T) 

can therefore be reduced to the problem of generating the function ee. 
Since la cannot be expressed as a rational fraction, the transfer func- 
tion for dead time can be simulated only approximately. Thus the remain- 
ing methods depend on an approximation of the transfer function gre by 
an approximate series or technique of network synthesis. 

4,2 Methods of Simulation 

4,21 Power Series 


The best known approximation to the function e is the Taylor 


series expansion 


ss X* Kies ye eX 
a asi 4,2 


This expansion is convergent for all values of x, but unfortunately 
the rate of convergence is small for x large. For this reason the 
Taylor series is not well suited to the generation of transportation 


lag for problems involving high frequencies or long time constants. 
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This can be seen by replacing s by jw in the expression for ae 


Moreover, in attempting to determine by longhand calculations, the trans- 
ient response of other than the simplest system possible involving trans- 
portation lag, the introduction of the Taylor series approximation for 


Ss - 
e f complicates the problem so Severely as to render it uSeless. 
4,22 The Pade Approximation 
In section 4.1 an important theorem of Laplace transforms was 


Stated, viz.: 
f[F(t-T) - Tle) es! repeated 4,1 


Thus the transfer function of transportation lag is 


_LLF-T)] ee oe ise 
Life] Ei 


Substituting s = jW to find the frequency response gives 


4.3 
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or the frequency transfer function of transportation lag is characteriz- 
ed by an amplitude ratio which is unity for all frequencies and a phase 
angle which varies linearly with frequency over the deSired range of in- 
terest, 

A popular method of approximating the transfer function "aie 
from this point of view is the Pade approximation. It stems from the 
fact that the quotient of any pair of complex-conjugate polynomials has 
an amplitude of unity and also possesses considerable phase shift. Ex- 


pressed in its basic form, the Pade approximation is 
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x LIM Fuv (x) 
G ae (u+V)—> a (one OF) 





4.6 


where VX Viv / x* 

v(x) = De vival) A — 
Fu.v (x) | (u+v) | é ce (Utv) (usv-1) 2! - 
(U+V) (U4V=1) © + of tt) VS 

Guwx) = f/- —4 ALLL 
i (Veuyl! * (V+u)(V+u-1) 2! 
-e by) ee ot 
(Vtu) (Veu-i)ees (vel) U! 
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The convergence for this series expansion is quite rapid and often values 
of u and v of 2 will give good accuracy for short time lags and low fre- 
quencies, The Pade approximant of enh for u = v = 2 is given in equa- 


tion 4.10 below. 
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In terms of the original problem, generating f£ (t-T) when f (t) 
is known, the equation for the second-order Pade approximant can be 


expressed in transfer-function notation as 


£(+-T) Be D*T*- GOT +!2 4.11 
f(t) D°T? +6DT +/2 


The justification of the validity of equation 4.11 is demon- 


strated as follows: 
L[ fT] = Foses™= [freer sereated 44 
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Substituting equation 4.11 in equation 4.1 gives 


LAT w STL este 
fle] ST*+6sT +2 


Rewriting equation 4.12 gives 
(?T'+6sT+2) Pf fle-7)] = (5:7  esT+l2) FL F(t) ae 
Performing the inverse Laplace transformation on equation 4.13 gives 
(D°T'+ GDT+12) f(t-T) = (D°T*-GDT +12) Fe) 4.14 


Comparison of equations 4.14 and 4.12 reveals that the equations are 
identical, This is, of course, a logical conclusion, as the Laplacian 
operator and the differential operator may be used interchangeably pro- 
viding the initial conditions in a System are identically zero. This is 
true in the system described above. 


The fourth order Pade approximant (u = v = 4) for = is 


ost _(ST) - 20(T)' + 18067) - 840 sT+ [680 
(sT)4 20(sT) + 180(ST)*% 84OST+ 1680 
The Pade approximants may be used for either longhand or com- 
puter solutions. For use in longhand solutions the approximant to be 
used is merely substituted for She It is worthy of note that the 
order of the characteristic equation is raised by the order of the 
approximant used to simulate a If the second order approximant is 
used, two additional roots are introduced into the characteristic equa- 
tion, and if the fourth order approximant is used, four additional roots 
are introduced into the characteristic equation. 


The technique of programming the Donner Analog Computer for 


handling the transfer function of the fourth order Pade approximant 
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is included in Appendix Section 9. 

The method of simulating transportation lag by the Pade ap- 
proximants yields what may be termed optimally flat phase characteris- 
tics at the origin but does not fit the phase characteristics as well 
over an extended frequency range. The major drawback to the use of the 
Pade approximants to simulate transportation lag is their response to a 
Step input or other input having high-frequency components, since the 
amplitude ratio is always unity but the phase is approximated only for 
lowWT products. The phase shift characteristics of certain of the Pade 
approximants are shown in Figure 4-2. As a result of the inability of 
these approximants to linearly phase shift high frequencies, the output 
of an approximant "rings'' badly for high-frequency inputs. This behavior 
is tllustrated in Figure 4-3 which depicts the response of the fourth 
order approximant to inputs of a step and ramp with unity feedback. It 
is noteworthy that the "ringing" is considerably less for the ramp in- 
put than for the step input. The reason, of course, is that the ramp in- 
put contains a lower percentage of high-frequency terms than does the step 
input. 

In the solution of problems concerning feedback control systems 
with transportation lag, advantage may be taken of the fact that the typi- 
cal system encountered is, in effect, a low pass filter. Block diagram 
algebra techniques are employed in the solution of these problems. Since 
the commutative law is applicable to block diagram algebra, the actual 
location of the transportation lag transfer function in the forward path 
is immaterial and it may be placed at the end of the path in the computer 
where the high frequency components have been attenuated to the greatest 


possible extent. Frequencies which have been attenuated 40 db. or more 
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are considered to be negligible and whether they have been linearly 

phase shifted is of no consequence. This tends to obviate the neces- 
sity of having to linearly phase shift a high WT product. Another 
consideration which makes unnecessary the linear phase shift of high WT 
products in the solution of the transient response of feedback control 
systems is that high w T products mean extremely large phase shifts and 
probable instability. In these cases the transient response is virtually 
meaningless. 

In the great majority of cases the fourth order Pade approxi- 
mant is more than adequate for the simulation of transportation lag. In 
cases where linear phase shift of higher wT products than those shown in 
figure 4-2 is desirable, higher order Pade approximants or more complex 
curve fitting techniques may be used to simulate aoe, The curve fitting 
techniques are discussed in the literature on analog computers but the 
need for such approaches is questionable. These techniques were devised 
because it was felt the Pade approximants of sixth order and higher fail- 
ed to phase shift low wT products accurately. In the course of this 
study the phase shift characteristics of the sixth and eighth order Pade 
approximants were calculated. The results are shown in figure 4-4, On 
the basis of these results it may be inferred that the order of the Pade 
approximation may be increased to increase the magnitude of the wT pro- 
duct that is linearly phase shifted without distorting the phase shift of 
low WT products. Because the Pade approximation may be expressed in a 
general form from which any order approximant may be developed, it is 
felt that the use of curve fitting techniques requires a needless expendi- 


ture of effort without any tangible benefit. 
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a. Effect of Time Delay on Linear System Analysis 


The addition of the time delay into the linear control system 
under investigation introduces several factors which tend to make the 
analysis of the system more difficult. Because the time delay may basical- 
ly be represented by an infinite series of terms, this creates an infinite 
number of roots. The effect of these complications on the standard graph- 
ical analysis techniques will now be investigated. In all cases, the 
system with a time delay of 1.0 seconds will be demonstrated. 

5.1 The Frequency Response 

The standard Bode diagram is shown in figure 5-1 when the time 
delay is included in the system. The difference between this plot and 
figure 3-1 is readily apparent. The phase margin has been reduced by a 
considerable amount. The gain margin is approximately 3.3 db which 
compares with the undelayed system which was stable at all gains. 

The time delay has introduced additional branches into the dia- 
gram. These branches appear due to the periodic nature of the time de- 
lay, and the increasing phase shift introduced at higher values of fre- 
quency. It is evident that the gain margin for each of the branches in- 
creases. This is analogous to the statement that the roots of the pri- 
mary branch would go unstable first if the gain were raised. 

The construction of this plot is simple and straightforward. 


For the system in question the open loop transfer function is 


~sT 
ec = — 5.1 
is S(0.1S+1) 


The gain vs logw curve is identical to the undelayed system 


curve since the time delay introduces no additional gain into the system, 
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The phase angle vs logw curve is drastically changed since 


the equation for the phase angle is 
i 
D = JO + ARC TAN 0,10 + W ID 552 


It is the presence of the additional term of the frequency times 
the time delay which introduces the additional phase shift and the extra 
branches. 

5.2 Root Locus 

The identical second order system with a transportation lag 
of 1.0 seconds also appears greatly changed when investigated by a 
root locus plot. 

Figure 3-2 showed the undelayed second order system which is 
stable for all gains. Figure 5-2 has the additional construction de- 
tail which we will require to obtain the root locus. In the same manner 
as the minus 180 degree phase shift line was obtained, lines of constant 
phase shift of 170, 160, 150, 140 degrees, and so forth, may be plotted. 
In addition, lines of constant phase shift caused by the one second time 
delay are also shown. In the non-dimensionalized plot shown, the time de- 
lay will produce a complete 360 degree phase shift every 0.628 radians. 

The resulting root locus is shown in figure 5-3. It is plot- 
ted by combining points of the basic system and the time delay which 
total to a 180 degree phase shift. It is evident that this system will 
go unstable with only a small increase in gain, and as a matter of fact, 
the gain at the stability limit just before the locus crosses the axis is 
approximately 1.48, the same value of gain margin obtained from the Bode 


diagram. 
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Figure 5-3. Root Locus of Second Order System With 
Transportation Lag of One Second. 
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An extremely interesting result of this plot is that upper 
branches of the root locus appear. Only two of these are shown although 
an infinite number exist. Here again they are a result of the 360 degree 
periodicity introduced by the time lag. It is also noted that at the 
same gain, the additional roots on the upper branches are attenuated more 
than the root on the primary branch, This substantiates the Bode diagram 
findings, which show that the roots on the primary branch would go unstable 
firse. 

A continuation of the graphical solution described above for 
other values of time delay will yield other root loci curves similar to 
the one in figure 5-3. 

The gain at the stability limit for each time delay was calculat- 
ed. The product of the time constant of the system and the gain, K T, , 
was plotted against the time delay on logarithmic cross-section paper. 
The result is a straight line as shown in figure 5-4. It was therefore 
possible to write a single equation which will describe the gain at the 
stability limit in terms of the time constant of the second order system 


and the time delay. 


k _ O.148 
‘Ca es 


Thus for second order systems if two of the above parameters are 


Tie 


known, the maximum value that the third parameter may have can readily be 


calculated. 


5.3 Transient Response 


Neither of the graphical methods discussed here will directly 


produce the transient response of any system. The Bode diagram will only 
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Figure 5-4. Non-Dimensionalized Gain at Stability Limit vs. 
Transportation Lag for Second Order System 
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yield values from which an approximate transient response may be infer- 
red. The root locus method provides the roots of the system on the s 
plane, but then the results must be transformed to the time domain and the 
equation of the output determined, The transient response must be obtain- 
ed by substituting values into this equation and obtaining the output at 
each specific time. Complicating this procedure is the fact that an in- 
finite number of roots exist on the root locus, and enough of these roots 
must be considered to obtain reasonable accuracy. Each time an addition- 
al set of these complex conjugate roots is considered, additional time and 
complexity are added to the above procedure. 

The transient response of this delayed second order system is 
shown in figure 5-5, It is here that the influence of the time delay on 
the second order system is most apparent. The original second order 
system was an overdamped system. It had no overshoot. The time delay has 
deteriorated the system to the point where the overshoot is approximately 
70 percent. 

It is interesting to note that a time delay of about 1.4 seconds 
would cause this system to be completely unstable. 

The effect of the time delay on the transient response there- 
fore is to increase overshoot and settling time. It can make a normally 
overdamped system highly oscillatory with a delay of only one second in- 


troduced, 
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6, Transient Response by Numerical Methods 


6.1 Discussion 

A mathematical technique now exists for the direct solution of 
the output of a linear continuous feedback system. An extension of this 
numerical method to the solution of systems with transportation lag is 
easily accomplished. 
6.2 The Z - Transformation 

A continuous time function f(t) is sampled by means of a sampl- 
ing switch S at equal sampling intervals of time T in figure 6-1. The 
output of the sampler f*(t) is a sequence of pulses which may be represent- 
ed for mathematical purposes as a series of impulses or delta functions 
whose areas are equal to the amplitude of f(t) at the respective sampl- 
ing ‘ad eons. Thus 


f(t) = Y $ (nT) §(t-nT) ae 


Ne ~¢o 


The Laplace transform of the sequence of impulses is given as 


tLriol= d (nt) 642 


+ 0d 
nso 


The Laplace transform of such pulsed functions, therefore, is in 


terms of a. Replacing this exponential term by a new variable Z, where 


Z =e nis accounts for the name Z-transform., One useful characteristic 
of the Z-transform is that the infinite series in equation 6.2 may usual- 
ly be expressed in closed form so that tables of Z-transforms may be de- 
veloped in the same manner as tables of ordinary Laplace transforms. 
One of the big advantages of the Z-transform is that it is possible 


to obtain the pulsed output of a linear system by expressing a pulse trans- 


fer function which relates the input and output Z-transforms. The pulse 
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Figure 6-1, Sampled Time Function 
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‘transfer functions can be expressed in closed form and tables of Z-trans- 
forms are available in most standard texts covering sampled-data systems. 
An abbreviated table is given in Table 6-1. 

6.3 Application to Linear Continuous Systems 

The application of the numerical method using Z-transform 
technique will be illustrated in a linear continuous feedback system as 
shown in figure 6-2(a). The system is mathematically converted into a 
sampled system by placing synchronized samplers at input and output and 
at point X in the feedback path in figure 6-2(b). The sampling frequency 
should be high enough to cause the error introduced by the sampled system 
to be small. By locating the sampler in the feedback loop and requiring 
that the transfer function G(s) be low pass in nature, then only the for- 
ward gain of the original system need be considered. Later examples will 
show the effect of the ratio of sampling frequency W, to the Bode cross- 
over frequency w, of the original system. 

Following the sampler S in the feedback loop is the element B(s), 
whose function is to reconstruct the original output of the forward path 
as faithfully as possible. This element reconstructs a polygonal approxi- 
mation to the actual function by converting each pulse into a triangular 
generating function as shown in figure 6-3. This requires the response of 
B to a short pulse or impulse of the same area to be a triangle as shown 
in figure 6-4(a). This response is not physically attainable since it 
requires a slope at negative time, but this has no effect on the purely 
mathematical solution. 

The combination of sampling switch operating with a period T 
and the triangle function generator B replace the direct connection 


which existed in the original system. The errors introduced into the 
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Figure 6-2, 3 Continuous Feedback System 
b) Sampled Model for Computation 
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Figure 6-4. ie Generating Triangle 
b) Ramps Which Produce Triangle 


37 





solution are a result of the inability of this combination to faithfully 
reproduce the feedback signal. 

The transfer function B(s) of the triangle function generator is 
obtained from the Laplace transform of the triangle function in figure 
6-4(a). This triangle is a result of the combination of three ramp func- 


tions as shown in figure 6-4(b). The Laplace transform of these three 








functions i8: 
a ea 
Bs) =e et eee 
or TS = 2 
R(sd= & | - ene 
ss ( . ) 


The relation between the Z-transform of the output and input in 
figure 6-2(b) is given by 


: Z = — OnG"(z) 6.4 
er) | + HBG*(2) 


and yields the pulse transfer functions for which tables are available. 

If the response of a linear system to an input F(s) is desired, 
first replace it with the system shown in figure 6-2(b). A proper choice 
of sampling frequency Ww, is made. The pulse transfer functions RG*(z) 
and HBG*(z) are computed, and combined according to the basic laws of 
algebra, yielding a form of a ratio of polynomials in powers of z. To 
obtain the output pulse sequence in the time domain, a simple process of 
long division is required. The division process will result in a power 
series of the form 


= i -3 es -~n 
O.(2) = Gas Guz) + Ouees: Qt +Oyé 4 och Cy 2 sae GTS 


The magnitude of the output in the time domain is just the coefficient 
of each of the terms in this series at a time corresponding to the 
power of z. Thus the output at the fourth sampling instant is ay, in the 


expansion, 
38 
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The long division process is continued until the desired 
number of points have been obtained. The effect of the sampling period 
To or of the sampling frequency UW, will be elaborated on in section 6.5. 
6.4 Application to Linear Continuous Systems with Transportation Lag. 

It has been shown in section 4.1 that the transfer function of a 
transportation lag or simple dead time lag is merely ey where T is the 
value Tp , or the time lag. 

If this transfer function is inserted into the linear system of 
figure 6-2(b) as G, (s) then the solution for the output will become 


2 OrGG, (2) 6.6 
Qe “(2) = loat- Sistalo- Gow ee) 


The value of this numerical method employing pulse transforms 
is now eaSily seen. The z-transform of the time delay ere F simply be- 
comes a) from the basic definition of the z-transform itself. Since the 
remaining terms of E)* (z) are also polynomials in powers of z, this 
does not add any more complexity to the basic long division solution of the 
transient response. 

It must be noted that the T introduced in the polynomial recon- 
Struction and the z-transforms in Table 6-1 is TS the sampling interval. 
This must be distinguished and kept separate from T 0 ,the time delay. 

One of the basic rules for the method is that the sampling time 
T. must be equal to or an integral sub-multiple of the time delay Ta 
It will be necessary to express ()*(z) only in terms of the sampling time 
Te in order that the value of the transient response will be obtained at 
each sampling interval. This method will be shown in the sample solu- 


tions in section 7. 
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6.5 Requirements for Reasonable Accuracy 

The ratio of Ww, to the crossover frequency WwW, of the open 
loop transfer function will determine the relative accuracy of the system. 
The comparison of sampling period T. with the expected settling time and/ 
or time delay in the system will determine whether the method will supply 
the required answer within a reasonable number of long division steps. 

Ideally, to obtain extremely good accuracy (within 1% of known 
answer) it is necessary to select a value of WwW, which is twice the fre- 
quency at which the open loop system is down 40 db on the Bode diagram. 
For the systems investigated here and explained in section 2, ™, is approxi- 
mately 1 rad/sec. For the system where G(s) = -, this would require WwW, 
= 200 rad/sec. Since Ds = a this would require a T. of approximately 
0.03 sec. With a normal settling time for this system of 4 seconds this 
would require a total of = or 133 long division steps to obtain the 
transient value at the end of 4 seconds. In order to reduce the computa- 
tional work and to remain within approximately an accuracy of 5%, vari- 
ous ratios of WwW, to W, were selected and results tabulated against re- 
sults of former traditional methods of solution. 

The problem of reducing computational work is further aggravated 
by the introduction of the time delay into the system. One basic re- 
quirement of this numerical method is that the sampling time ia must be 
equal to or an integral sub-multiple of the time delay T, . For systems 
with a long settling time, for examp le 5 seconds, and a short time delay 
of 0.10 seconds, it would require 50 long division steps to produce the 
value of the transient response at 5 seconds. This is only if it would 


be possible to select T. = T In the event that the value of TS must 


0 e 


be chosen smaller to satisfy the requirement of the ratio W, to w, for 
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accuracy, then it would require 50n long division steps, where 





It will be demonstrated in the next section that this method 
will provide reasonable engineering accuracy (approximately 5%) with a 
small investment of time. This can be accomplished with a ratio of WwW, 
to Ww, as low as possible, and a T. which is the lowest sub-multiple of 
T D as long as it is not greater than = : 

This method can give extremely good accuracy at the expense of 
an extremely long and repetitious series of calculations. Its advantage 
lies in the fact that where the ratios mentioned above are favorable, it 
should be possible to compute the transient response of a system in a 
short period of time with no more powerful tools than algebra and a slide 
rule or desk calculator. 

The next section will demonstrate the accuracy and simplicity 
of the method in solving simple feedback systems. Computational time 


and labor will be compared with other methods currently available to the 


engineer. 
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ye Comparison of Numerical Method with Other Methods of Solution 


7.1 Numerical Method Solution 

The numerical method solution to the second order system des- 
cribed in section 2 will now be demonstrated in detail. The system to 
be solved is shown in figure 6-2(b). The input is a unit step. The time 
delay will be one second. The sampling time interval will also be one 
second. The table of results in section 7.5 will also show the solution 
results for other sampling intervals. 


The output ©. (s) in the Laplace notation is 


/O Sse 


_ S*(S+10) 
T SC S+10) 


where 


er 


Bes ol — ee 


sl) 
oa ) 


deme 


Combining 10 ef 


2 S* (5+/0) 
O.(s) =  Jotiaeeer Pea 


[ + S57 Sey 


It is important to note that the T in the time delay represents 
Tp , the delay time, whereas the T in the expression for B(s) and also in 
the Z transforms represents To the sampling time. In the algebraic work 
in combining the expressions it will be necessary to express Tp in terms 
of a and obtain the output equation in terms of T. only. In this ex- 
ample, since Tp equals 7. equals one second, it will not be necessary to 
account for this difference. 

By expressing numerator and denominator separately in partial 
fraction expansions, and then transforming each separately by means of 


the transforms in table I of section 6 we obtain 
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O.(2) = ys 
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200(i-z") 


Combining numerator and denominator results in 


Q.(2) = 1$0z*+ 2027 oe 
200 - 318 27'+2342°7-1/427>-22 
By carrying out the long division indicated in equation 7.5 


we obtain a power series in Z 


Bees 0.9274 1532 + 1 Beer Og) ee a 


We may directly transform this expression into the time domain. 
It is not really necessary since the only effect on the output will be 
to express the series in terms of the sampling periods. [In either case 
the coefficients are the same and are the value of the transient response 


at each sampling instant. 
Q.(t) = 0.9(t-2T) + 1. S31(t-37) + 1.385 (t-4T) + °° jeg 


The introduction of a different time delay or sampling time 
into the problem merely changes the algebraic relationship between Tp 
and T. and once this accounting has been accomplished, the algiaitte and 
final solution work is similar to that shown here. No additional diffi- 
culty is introduced into the problem, 

7.2 Power Series 
7.21 First Order System 

The solution of the transient response of a first order system 

with transportation lag by the employment of a power series expansion ia 
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quite straight forward, although tedious and time consuming. The 
System under consideration is that shown in figure 7-1l(a) where T is 
equal to one second. 
The closed loop transfer function may be written 
cn 
fees — _ 7.8 
Op {+ o° 


For & unit step input 
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The quantity ( l+e ) may be expanded in a power series 


with the following result. 


~TS -Ts -21s ~3TS ~4TS 
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Taking the inverse Laplace transform of each term yields 
O<(E) aa eee i) ~ ~ u(t-2) fe 2) ete JSP 3)" 


7.14 
— jr (t-4) {eo + w(t-5) te" = 
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Figure 7-l(a). First Order Feedback Control System 
With Transportation Lag. 
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Pigure 7-1(b). Second Order Feedback Control Systen 
with Transportation Lag. 
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The transient response may now be determined by substituting 
values of t into equation 7.14 . The magnitude of the output at any 
time t may be determined as exactly as desired by evaluating all terms 
and by evaluating each term to a sufficient number of decimal places. 
The amount of time required to obtain a transient solution by this method 
increases directly as the time interval between sampling instants is de- 
creased, It also increases as the settling time increases since increased 
settling time increases the number of terms which must be evaluated at 
higher values of t, The sample problems solved by this method required 
about twelve hours to complete. 
7.22 Second Order System 

While the use of a power series expansion is a straight forward 
method of determining the transient response of a first order system with 
transportation lag, the method requires far too much work to be considered 
practical for higher order systems. The following example will demons- 
trate why the method fails-as a practical working tool. The system under 
consideration is that shown in figure 7-1(b) where T is equal to one 
second. 


The closed loop transfer function may be written 


eee 
Oe a SCO1S +1) 
Oe 7 ens 7.15 
a S(0.1 $+1) 
For a unit step input 
QO. = species 
. S 
Then est 7 
$* (0.1 S41) ee poi : 
— — ar eee 7.16 
ee $*(0./s+1) S(Oj/S+1) 
S (0.1 S41) 
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may be expanded in a 


est =| 
50.15 +!) ) 
power series with the following results. 


The quantity [1+ 


-sT -sT -2sT 
=  @ ¥.* a ce + e 
S*(O.1541) SCO.1S +t) S*(0.15+1)* 
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a e 
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SAlalscol.) S"(O.4S+1)* 


At this point it may readily be seen that the task of determin- 
ing the inverse Laplace transform of each term is momentus and the work 
involved increases tremendously with each increase in the order of the 
term. Thus it is maintained that the use of a power series is impracti- 
cal in determining the transient response of systems of second order or 
higher which are complicated by transportation lag. 

7.3 Pade Approximation (Arithmetic) 

The transient response will now be determined for the system 
of section 7,1 by simulating the Laplace shift operator with the second 
order Pade approximant, The major steps in the solution will be outlined 
but the tedious arithmetical calculations will of necessity not be repeated 
in their entirety. 

The system in block diagram form is shown in Figure 7-1(b). 
Since the duration of the transportation lag is one second the Laplace 


shift operator becomes e°, As indicated in section 4,22 


rc Se ee 





7.18 
. §*+6S +12 
Substituting the Pade approximant in for the Laplace shift 
operator and manipulating by block diagram algebra yields the closed 
loop transfer function for the system, 
6 lo(S*%~6§+12) 
. = S(S+10) (57465412 — 
Or | + 10(€S*~GS +12) ‘ 
S Cor os 46s 4712) 
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For a step input, OQ, = ge 


S 
i T 1OS?>—6OSe+ 120 ) 
— ree 2 
Oc S ne ee. ae 


The transient response of this system to a step input is con- 
tained in Laplace notation above. The problem of reverting to the time 
domain remains, 

The first step is the determination of the four roots of 
st + 16S> oe 828° + 60S + 120. Since these roots are two complex conju- 
gate pairs whose real and imaginary parts are not necessarily integral 
the Lin method was selected for their solution. The first trial divisor is 
828° + 60S + 120 which equals 57 + .7338 + 1.46. Five iterations were 
necessary to yield factors of ( 5? + .4638 + 1.65 ) ( 5° + 15.548 \+ PSiedu 
Application of the quadratic formula further reduced the factors to ( S + 
0232 + j1.264 ) ( S + 7.77 + j3.59 ). 

Next, the method of residues was employed to determine the co- 


efficient of the inverse transform of each of the five poles. 


Thus 


pee (277443 s59t 


Q. (ec) = [+ (G32 £4.45) + (</8/ 24.232) € ath 


This further reduces to 


-,232U 
Oc (t) = i =e ee i (2 cosh2e4 e+ 3.46 sin /.264 €) 


= Fie 
~ AH” (2 cos 3.59E + 2.63 sin 3.59€) he? 


The full transient response of the System may now be plotted by 


evaluating Oe (t) at various values of t. 
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The time required for the solution of the transient response 
of the system by this methoc is between seven and eight hours, It is 
worthy of note that the roots of the characteristic equation determined 
by this method are reasonably close to the values obtained from the more 
accurate but more lengthy root locus solution. Arithmetic solutions 
making use of the Pade approximation were limited to the use of the 
second and third order approximants because the time required to solve 
a probiem when higher order approximants are used becomes unreasonably 
long. 

7.4 Pade Approximation (Computer) 

In section 4.22 it was pointed out that the fourth order Pade 
approximant was more than adequate for handling the vast majority of prob- 
lems involving transportation lag. Accordingly, all problems solved by 
arithmetical techniques were also solved on the Donner Analog Computer 
using the fourth order Pade approximant to simulate epee. The solution 
to the problem by this method was established as the control against 
which the arithmetic solutions were compared, The fourth order approxi- 
mant was chosen because it occupied only six amplifiers on the Donner 
Computer thus leaving four amplifiers free on which the remainder of the 
system could be programmed. The computer program for the fourth order 
Pade approximant is contained in the Appendix. This program was specifi- 
cally written for a range of values of transportation lag between 0.1 
seconds to 1.0 seconds. In each amplifier only one resistor R. had to 
be changed to change the value of the transportation lag between the limits 
Stated above. For transportation lags outside of these limits the program 
should be modified so that Ry does not exceed the bounds of 0.1 to 1.0 


megohm. 
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7.5 Table of Results 

This table of results is an intent to show the relative ac- 
curacy of the various methods, and to compare the solution time in- 
volved to achieve this accuracy. The first order system described in 
section 2 was used to establish the basic comparisons for accuracy, and 
to demonstrate the feasibility of the numerical method. Results indi- 
cate that the numerical method approach provides results of comparable 
accuracy to other methods with substantial savings in computational time. 
Results also indicate that increased accuracy may be obtained with the 
numerical methods solution by increasing the sampling frequency. It 
must be pointed out, however, that the time to arrive at the solution 
increases directly in proportion to the increase in the sampling fre- 
quency. 

Table 7-1 contains the basic information from which general 
conclusions about relative accuracy will be based. The results indicate 
the value of the transient response for a first order system with a one 
second time delay as obtained by various methods, Results are shown at 
one second intervals. The power series solution is the most accurate 
available and will be taken as the standard by which the other methods 
will be compared. 

Table 7-2 contains similar information for the first order sys- 
tem with a 0.5 second delay. The numerical method solution for a ratio 
of Ws to Up ratio of 12.56 is compared because it achieves the stated 
requirements of the study of less than 5% deviation. It will be the ratio 
used in higher order systems. 

Tables 7-3 and 7-4 display results for the second and third order 


Systems respectively for a 0.5 second delay. The basis for comparison in 
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these systems was the computer solution employing the fourth order Pade 
approximant. The deviations here from the standard are of the same 
orders of magnitude as experienced in the first order system. In these 
results, a plus deviation indicates a result higher than the standard, 
and a minus deviation a result lower. The numerical method solution 

when compared with the standard characteristically maintains the same 
sign at the same relative point of the transient response. That is, the 
transient response as determined by the numerical method coincides with 
or leads the standard, and would therefore tend to show a slightly short- 
er rise time and time to peak overshoot. Its accuracy in determining the 
peak overshoot and settling time was well within the engineering accuracy 
prescribed for this investigation. 

It must be understood that the solutions by the numerical methods 
technique yield "exact'’ values only at the sampling instants. Interpola- 
tion may be done bo determine information at times between these sampling 
instants with assurance of reasonable accuracy as the magnitude and alter- 
nating signs of the deviations of the tabulated results from the standard 
indicate. 

All solutions, both computer and analytic, employing the Pade 
approximation exhibit large deviations at the first few tabulated times. 
These are due to the "ringing" discussed in section 4.22 and illustrated 
in figure 4-3. This is an unavoidable disadvantage of the Pade approxi- 
mation of . 

It cannot be over emphasized that the time required to obtain 
the transient solution by the numerical method was in most cases from 


one-fifth to one-tenth the time required by previously used analytical 
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and graphical methods. Its value also lies in the fact that the pro- 
cedure is simple and easy to follow. Once the engineer or student 
becomes familiar with the procedures, it is possible to recognize and 
correct computational errors before they are compounded in following 
steps. This is simply not possible in the other analytical methods 
until the final step of the solution. Similarly, a change in the 
transportation lag or the sampling interval is easily handled, and 

a great deal of the algebraic manipulations need be done only once 


for each basic system, 
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COMPARISON OF TRANSIENT RESPONSE AS OBTAINED BY 


TABIE 7-1 


POWER SERIES WITH OTHER ARITHMETIC METHODS 


FIRST ORDER SYSTEM - TRANSPORTATION LAG EQUAL TO ONE SECOND 


Power 
Series 


1.00 
1.50 
1.167 
f91 
-8416 
1.057 
1.100 
1.004 
948 
9798 


1.021 


0499 
848 
1.055 
1.098 
1,005 
950 
981 


1.019 


+1.3 
-1.0 
-0.8 


+0.5 


+0.2 
+0.1 


-0,2 


1.00 
1.50 
1.250 
.875 
8125 
. 9688 
1.0782 
1.0548 
9884 
9669 


9894 


0 
+7.1 
+10.6 


=o 


+4.9 


+4,2 
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Num. 
Meth. 
0 
0 
1.00 
1,50 
1.1875 
~8125 
8820 
1.0331 
1.0993 
1.0205 
© 9564 
~9719 


1.0108 


Ws 


+1.8 


+2.7 


ee 


-2.3 


-0.1 


+1.6 


+0.9 


-0.8 


-1.0 


Num. 
Meth, 


1.00 
1.50 
1.1788 
- 7969 
8392 
1.0509 
1.1022 
1.0087 
sa943 
9916 


1.0329 


“i 
Dev. 


+0.41 
+0.75 
-0.28 
=0736 
+0.2 

40.37 
+0.66 
+1.20 


41.15 
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TABLE 7-2 
COMPARISON OF TRANSIENT RESPONSE AS OBTAINED 
BY POWER SERIES WITH NUMERICAL METHOD RESULT 


FIRST ORDER SYSTEM-TRANSPORTATION LAG EQUAL TO 0.5 SECOND 
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Power Numerical we 
Metho 
Series Dev. 
icon 
ioe = 12.56 
0 ¢) 0 
0 0 0 
0.5 0.5 0 
.875 2875 0 
1.028 1.031 40.29 
1.039 1.055 +1.06 
1.021 1.033 +1.18 
1.006 1.011 +0.5 
999 1.000 +0.1 
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TABLE 7-3 
COMPARISON OF TRANSIENT RESPONSE FROM COMPUTER SOLUTION 
WITH NUMERICAL METHOD AND PADE APPROXIMANT (ARITHMETIC) 


SECOND ORDER SYSTEM - TRANSPORTATION LAG EQUAL TO 0.5 SECOND 
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Computer Numerical uF Pade Te 
Method Approx. Dev. 
_— Dev (Arith) 

0 0 - sO. - 

0 0 - 0018 - 
wal aD +12.8 ~404 +30.4 
77 803 + 4.3 846 + 9.9 
1.04 1,078 + 3.7 1.073 + 3.2 
Las 1.153 + 3.4 1a - 0.5 
1.100 1.110 + 0.9 1.066 - 1 
1.050 1.037 - 2 1.017 = 3et 
1.015 984 ~ Sal 992 = 2S 





TABLE 7-4 
COMPARISON OF TRANSIENT RESPONSE FROM COMPUTER SOLUTION 
WITH NUMERICAL METHOD AND PADE APPROXIMANT (ARITHMETIC) 


THIRD ORDER SYSTEM - TRANSPORTATION LAG EQUAL TO 0.5 SECOND 
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Time Numerical Pad y 
Computer Method to ’ 
Sec, Dev Approx. Dev 
H-12,56 (Arith) 
0 0 0 = ® 146 i 
0.5 0 0 as ~0203 = 
1.0 26 aie +34.6 » 359 +38.0 
1.5 hed ~803 +10.0 821 +12.5 
2.0 1.050 1,078 + 2,7 1.09 + 3.8 
2.5 1,150 L, B53 + 0,3 1.150 0 
3.0 beob15 1.110 - 0.5 1.099 - 1.4 
a, 5 1.060 1.037 - 2,2 1.029 - 2.9 
4.0 1.005 984 - 2.1 987 - 1.8 
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8, Conclusions and Recommendations for Further Investigation 


The investigation into the effect of transportation lag on 
feedback control systems reported on in this thesis has demonstrated 
four ideas of practical significance. It has also opened up a number 
of avenues for further investigation. 

The primary conclusion which the investigation supported is 
the applicability of the numerical method technique devised by RAGAZZINI 
and BERGEN to the problem of determining the transient response of linear 
feedback control systems which are complicated by transportation lag. Of 
all the mathematical and graphical tools available to the engineer it 
offers the optimum combination of speed, accuracy, and simplicity. It 
may be used by students whose background includes courses only as advanc- 
ed as linear feedback control systems, provided the Z-transform is included. 

The deviation of the transient response as determined by the 
numerical method from that determined by the techniques currently in use 
may be held to less than three percent. Meanwhile the time required to 
determine the transient response has been reduced by more than seventy- 
five percent. Thus the decision of whether compensation is required or 
not may be made much more rapidly than before. The form the compensation 
will take will, of course, vary from problem to problem and from designer 
to designer. 

The suitability of the numerical method to the particular pro- 
blem of transportation lag suggests an area of further investigation. In 
determining the transient response the system has been mathematically con- 
verted to a sampled-data system, and is described by Z-transforms. In 


order to determine whether the system will meet specified performance 
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criteria it is necessary to calculate the transient response. If the 
System is not satisfactory, compensation techniques applicable to con- 
tinuous feedback systems would be employed and the transient response of 
the compensated system recalculated. However, since the system is al- 
ready a mathematical sampled-data system, it would appear that the suit- 
ability of sampled-data system performance criteria and compensation 
techniques to the problem is worthy of investigation. Should these 
techniques prove successful the reverting back and forth from the -Laplace 
transform to the Z-transform for compensation could be avoided and a 
greater saving in time and tedium effected. 

The second noteworthy conclusion of this investigation is 
that the Bode diagram not only offers the fastest tool for the deter- 
mination of system stability but it also will indicate which set of roots 
will go unstable first. This information is made available by no other 
technique except the root locus and a conservative estimate would place 
the time differential using the root locus at a full order of magnitude 
higher than the Bode, 

The criterion of dominance of a set of roots in a higher order 
system provides the basis for the extension of the performance character- 
istics of a second order system to these higher order systems. From the 
information in the Bode and root locus graphical plots, it may be practical 
to construct a set of second order curves with transportation lag which 
might have applicability to the analysis of higher order systems. 

A third significant conclusion is that the Pade approximation 
for the Laplace shift operator is considerably more linear in its phase 
shift characteristics than the recent analog computer literature would 
indicate. An inspection of figure 4-4 will reveal that the Pade approxi- 


mation provides a completely linear phase shift until that point where it 
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deviates from the ideal, Thus, the simulation of time delay on the analog 
computer requires merely the choice of a high enough order Pade approxi- 
mant which will linearly phase shift the frequencies desired. It would 
appear that the curve fitting techniques currently in existence are not 
justified due to their increased complexity and lack of a simple repeti- 
tive scheme for programming, 

A fourth interesting result of this study is the information on 
stability limits for the second order system. When plotted on logarithmic 
cross-section paper as shown in Figure 5-4, the non-dimensionalized gain 
vs. the time delay produced a straight line curve from which it was possible 
to write an equation describing the relationship between gain at stability 
limit and transportation lag for any second order system. 

This result is obtainable from root locus or frequency response plots 
and is easily presented because of the ease in non-dimensionalizing a 
second order system. Here again, a further extension of this informa- 
tion would be possible for higher order systems, The method of display- 
ing the information would depend on the way in which the additional time 


constants of the higher order systems were handled. 
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9. Appendix - Computer Program for Fourth Order Pade Approximant 


The Pade approximation of e* is discussed in Chapter 4, An 
analog computer program for the fourth order approximant which was used 
in the study is developed here. 


The Pade approximation states that 


eT Wien Gy e 
where 
=a) = | + aah aaa 
(u+v)(Utv-l)++- (UFl) VI a 
and 
Gm) = | - Xo. U(u-l) Xe, 
| (v+u) 1! (Vtu) (vtu-t) 2! 
eG! (in eae me 
Substituting u = v = 4, and x = -sT yields 


e . (ST)*- 20(sT)*+ /2O(ST)*- 840(ST) + I6FO s 
(ST) "+ 20(sT)*+ 1 90 (ST)*+4 840(ST) +1680 


20 /8O SHO {GSO 
ee. =) 6 on (ee orn) 
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As shown in Section 4.22 


Pw LEC} 9.6 
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Since the differential operator and Laplace operator are inter- 
changeable when initial conditions are zero, equations 9.5 and 9.6 may be 


combined as follows 


j~ 29 + 180. _ Sto 1680 





f(tt) 2 __/ opr pert ps t 2 
t 20... “80 * Sao [G80 
F(t) f+ 37+ pret pep t pep 


The final form of equation 9.7 from which the analog computer program 


shown in figure 9-1 was developed is 


(e-T) = $e) - (29°F + 22 4¢cer] 
180 = 
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TABLE OF COMPUTER COMPONENTS 


Potentiometer Settings Component Values 
ay = Ro, 1 Megohm 
: ay = RoC, RR Ry 
ay = RAC, R, Re 
ay = 0,2 Co Re Re 
a 
5 = Ree» 5 Ri 
268 laecs Eee 
ay = 0.467 Cc, Ree 
a. = Rec, 
Ag = Ry, 2 Megohm 
a1 =z 0,9 C), Res 
Sage) yee 
R 
412 * 1 1 Microfarad 
£5 
2. 
a Se ele C C 
is Res 1 2 
R 
£6 


The time is scaled such that computer time is equal 


to ten times problem time. 
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